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Abstract. Several doped semiconductors, in contrast to heavily-doped silicon and 
germanium, host extremely mobile carriers, which give rise to quantum oscillations 
detectable in relatively low magnetic fields. The small Fermi energy in these dilute 
metals quantifies the depth of the Fermi sea. When the carrier density exceeds a 
threshold, accessible thanks to the long Bohr radius of the parent insulator, the local 
seafloor is carved by distant dopants. In such conditions, with a random distribution 
of dopants, the probability of finding an island or a trench depends on on the effective 
Bohr radius, a^ and the carrier density, n. This picture yields an expression for electron 
mobility with a random distribution of dopants: oc(a n in reasonable 

agreement with the magnitude and concentration dependence of the low-temperature 
mobility in three dilute metals whose insulating parents are a wide-gap (SrTiOa), a 
narrow-gap (PbTe) and a “zero”-gap (TlBiSSe) semiconductor. 


1. Introduction 

In an outreach paper written in the heyday of elemental Fermiology, Mackintosh 
remarked that while “few people would dehne a metal as a solid with a Fermi surface, 
this may nevertheless be the most meaningful dehnition of a metal.” [1]. In opposition 
to a metal, a semiconductor can be dehned as a solid deprived of a Fermi surface. 

Doping a semiconductor [2], can eventually turn it to to metal. This metal-insulator 
transition [3] has been a central issue of the condensed-matter physics during the second 
half of the twentieth century. The most-explored system has been silicon doped with 
phosphorus, its immediate neighbor to the right in the periodical table. Experiment 
has identihed a sharp transition and has accurately measured the critical concentration 
of dopants for emergence of metallicityjl]. Above this threshold concentration, silicon 
becomes metallic: It conducts electricity at zero-temperature. This is an alternative 
dehnition of a metal[5] and more inclusive than the one mentioned above. Historically, 
the debate on metal-insulator transition have been focused on “weak metallicity”, which 
only requires the presence of mobile electrons at zero temperature and not “strong 
metallicity”, which implies the existence of a well-dehned Fermi surface. In a variety of 
insulators, the critical concentration of dopants for metal-insulator transition has been 
identihed by experiment |6] and the results match the expectations of the Mott criterion 
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Figure 1. Quantum oscillations of the Nernst coefficient (a) and the temperature 
dependence of the amplitude of the oscillations (b) in SrTiOa-i [adapted from ref. [12])]. 
They attest to the existence of a well-defined Fermi surface and quantify the amplitude 
of the Fermi energy, c) The results imply that the chemical potential is ImeV above 
the bottom of the conduction band separated from the valence band by a gap of 3 eV. 
d) The seafloor is this shallow Fermi sea is rough and the low-temperature mean-free- 
path of electrons is set by the average distance between islands or trenches, which are 
potential scattering centers. 


for metal-insulator transition [Hi, which identihes the effective Bohr radius of the host 
insulator as the key parameter in a given material. 

According to the Pauli exclusion principle, electrons of a solid which share the 
same discrete quantum numbers (such as spin and valley) cannot have the same wave- 
vector. They will fill k-states in a Fermi-Dirac distribution with an edge separating 
the hlled and empty states. The sharpness of this edge is set by the inverse of 
the mean-free-path, which defines the size of pixel in the k space. A sharp Fermi 
surface has an experimental signature other than a finite zero-temperature electric 
conductance. Landau quantification truncates a Fermi surface to a number of concentric 
tubes leading to quantum oscillations in various physical properties [9]. In many doped 
semiconductors, such oscillations have been reported decades ago. The list includes 
n-doped strontium titanate[Tni IHl [I2l US] [H], doped IV-VI semiconductors such as 
PBS [ig, PbTedSl HSl CZ] and SnTe[IHl US]- This is also the case of heavily-doped 
germanium [201 [2T| 122] . A recent study of Nernst effect in SrTiOs found quantum 
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SYSTEM 

e/eo 

m* /me 

a^ (nm) 

Si 

12.5 

0.45 

1.5 

Ge 

16 

0.24 

3.5 

SrTiOs 

20000 

1.8 

600 

PbTe 

1000 

0.07 

800 


Table 1. Normalised electric permittivity, e, the average effective, to*, and the effective 
Bohr radius (a’Jj = ), in four different semiconductors. 

oscillations of a very large magnitude with a single and well-defined frequency (See 
Fig.l). This result implies that this wide-gap semiconductor, after loosing one out of 
its 10^ oxygen atoms, becomes a metal in the strongest sense of the word: a solid with 
a Fermi surface. 

From this perspective, one unexpected hnding of the recent research on topological 
insulators is instructive. The quest to find ultra-mobile electrons on the surface 
protected from backscattering mostly documented bulk mobile electrons giving rise 
to easily-detectable quantum oscillations [23l IMl ESI ES]. It is now clear that 
topological insulators are often on the metallic-side of the metal-insulator transition 
as a consequence of uncontrolled doping and a long Bohr radius in the bulk. The irony 
was not pursued however. Bulk electrons were often treated as a nuisance, impeding 
topological “insulators” to qualify as true insulators. Their proper mobility, with notable 
exceptions |27], was not put under scrutiny. 

The subject of the present paper is the intrinsic limit to carrier mobility in a dilute 
metal made by the introduction of a random distribution of dopants in a semiconductor 
with sufficiently long Bohr radius. A pattern is visible in the available experimental 
data. In these systems, at identical carrier density, mobility of electrons is much higher 
than in metallic silicon or metallic germanium. Moreover, the variation with carrier 
concentration is quite distinct. The main argument of this paper is that both the 
magnitude of the mobility and its density dependence can be understood by conceiving 
the simplest model for the probability of hnding a scattering center in a shallow Fermi 
sea. The key parameter remains the Bohr radius of the host insulator, which by setting 
the Thomas-Fermi screening length, shapes the local seahoor of the Fermi sea. 

2. Carrier mobility in metallic semiconductors 

Fig. 2 Compares the Hall mobility at liquid-Helium temperature of of four different 
semiconductors as a function of carrier concentration between 5 x and 

10^^cm“^. The reported data for phosphorus-doped silicon [28], arsenic-doped [29], 
antimony-doped [30] and gallium-doped]^ germanium are presented together with those 
for p-dop edm m and n-doped|32] PbTe, as well as for oxygen-dehcient SrTiOa 
[H E31 El ES]. The hgure reveals several remarkable features. 

The hrst is the magnitude of the low-temperature mobility. The room-temperature 
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Figure 2. Mobility as a function of carrier density in a number of doped 
semiconductors. Neither doped silicon, nor doped germanium show a strong variation 
in mobility, when the carrier density is above the threshold of metal-insulator 
transition. In both cases, the metal-insulator transition and the onset of localisation 
(corresponding to ~ 1 occurs almost simultaneously. On the other hand, in PbTe 
and SrTiOa-^, which have both a long Bohr radius, mobility enhances with decreasing 
carrier concentration and localization is easily avoided. Solid horizontal lines mark the 
required mobility to see detect quantum oscillations at a given magnetic field. 


mobility in both SrTiOs or PbTe is disappointingly low, closing the door to numerous 
applications. But as electrons are cooled down the resistivity decreases by several orders 
of magnitude. The amplitude of inelastic scattering, a remarkable subject by itself|HB] 
will not be discussed here. The focus of the present paper is the upper bound to mobility 
set by elastic scattering, which as one can see in the hgure, is very different among the 
four semiconductors. At a carrier density of the low-temperature mobility 

in PbTe is orders of magnitude larger than in SrTiOs, which in its turn exceeds by 
far the mobility in germanium. At this carrier density, silicon is an insulator and its 
low-temperature mobility has dropped to zero. 

There is a second striking feature in the hgure. While metallic silicon and 
germanium present an almost hat mobility in the concentration window just above their 
metal-insulator transition, in both PbTe and SrTiOa, the mobility clearly varies with 
carrier concentration. Remarkably, this allows them to keep a large ]^p^ in the range 
under consideration and to avoid localisation. On the other hand, the hat mobility in 
silicon and germanium makes them vulnerable to localization as ki? is reduced. In both 
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cases, the onset of metal-insulator transition (which matches the Mott criterion) occurs 
close to where the condition for Anderson localisation ~ 1) is satisfied. As Mott 
himself remarked in the preface of the 1990 edition of his book|3], while metal-insulator 
transition in Si;P occurs at the critical density corresponding to the Mott criterion, it 
is probably an Anderson transition driven by disorder. 

Finally, the figure instructs us on the amplitude of magnetic field required to 
observe quantum oscillations. Roughly one needs fiB > 1 to detect them. It is 
not surprising that there is no trace of a report on quantum oscillations in heavily- 
doped silicon, given its mobility. The large electron mobility in dilute PbTe is the 
reason that an electromagnet providing a field as small was sufficient for high-resolution 
fermiology[l6l [T7]. In the case of STiOa, quantum oscillations become more pronounced 
as carrier concentration is reduced and the insulating state is approached |12j. This 
would only look like a paradox if one forgets the density dependence of carrier mobility. 

What distinguishes SrTiOa and PbTe from silicon and germanium is their 
exceptionally long Bohr radius. The relevant parameters for each system are summarized 
in Table 1. A Bohr radius approaching micron explains why both systems remain strong 
metals at a carrier density of when silicon and germanium have already 

become (or are still) insulators. The purpose of this paper is to consider the consequences 
of a Bohr radius deep in the metallic state. 


3. Three relevant length scales 


Decades ago, Mott began to draw a picture of how a critical concentration of dopants 
transforms an insulator to a metal[71 18]. In his argument, three length scales play a key 
role. 

The first one is the effective Bohr radius of the host insulator. The Bohr radius was 
introduced in the early days of the quantum mechanics as a length scale connecting four 
different fundamental constants, the Planck constant,/!, the mass, m^, and the charge, 
e, of electron and the electric permittivity of vacuum, cq: 

o-B = = 0.053nm (1) 

The Bohr radius quantifies the characteristic length scale of the Coulomb interaction 
in a hydrogen atom. Now, in a crystal, two out of the four parameters can differ from 
their magnitude in vacuum. The effective mass of quasi-particles depends on their 
energy dispersion. The electric permittivity (the dielectric constant) is set by the way 
the medium responds to an applied electric field. The effective Bohr radius reffects these 
features: 


e me 
Qb = -dB 
eo m* 


( 2 ) 


Thus, a large dielectric coefficient and a light effective mass elongate the Bohr radius 
and enhance the spatial scope of the electric field generated by an ionized dopant. 
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Figure 3. a) and b) Illustration of the Mott criterion for the metal-insulator transition. 

c) When the Bohr radius becomes much longer than the interdopant distance, Bohr 
spheres interpenetrate. Scattering is governed by statistically unavoidable regions 
(marked by arrows) in which the local concentration is well below the global average. 

d) The two relevant length scales in the metallic side of the transition are the screening 
length, vtf , which determine the spatial decay of the potential energy around a dopant 
and the interdopant distance, _ g) xhe local energy landscape is set by the sum of 
all contributions from all potential wells, f) Poisson distribution becomes wider as A, 
the expectancy, increases. For A > 10, it can be approximated as a normal distribution. 
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The second relevant length scale is the Thomas-Fermi screening length, ttf, which 
sets the length scale for the spatial decay of a screened Coulomb (Yukawa) potential, 
which is expressed as : 


V{r) = e 


exp{—r /rTp) 

4ner 


(3) 


Thomas-Fermi screening describes the response of a Fermi-Dirac distribution to a 
perturbation in charge distribution. ZimanPT] showed that using the Poisson equation 
for a perturbative electric potential, one can the following expression for ttf (expressed 
here in S. I. units): 


2 ^ ^0 _ ^ 

em{eF) 4kF 

Roughly speaking, while a^ quantihes the spatial extension of a Coulomb potential 
well in an insulating medium, i^f is a measure of how spatially fast the well is screened 
in a metal. The Mott’s criterion for metal-insulator transition corresponds to: 


rTF — Cl*B 


(5) 


When Ttf becomes longer than the screening of Coulomb interaction in the 
metal extends over a distance longer than the size of the well in the insulting medium. 
Electron-hole pairs become unstable and electrons trapped in adjacent potential wells 
connect and form a conduction band. Ashcroft noticed that this is the critical value for 
the emergence of a bound level in a Thomas-Fermi held and this suffices to make the 
Fermi sea unstable |38]. 

Mott expressed his criterion in terms of the interelectron distance, which is n“^/^ 
in isotropic medium with a density of n. This is our third length scale, which is linked 
to the Fermi wave-vector in a Fermi-Dirac distribution : 


n = 


7,3 

Stt^ 


( 6 ) 


One can use Eq. 2, 4 and 6 to rewrite Eq. 5 and derive the usual expression of the 
Mott criterion: 


1 


a^n^/3 = -(-)^/3 _ q_253 


4^3 


(7) 


The Mott criterion has been derived in different ways. Using a polarization- 
potential method, Bhatt[39l SO] showed that in a cubic lattice of hydrogen atoms, the 
gap between upper and lower Hubbard bands closes when the lattice constant is reduced 
to approximately four times the Bohr radius. In absence of disorder, the transition is 
expected to be hrst order. 

Experimental data on a variety of systems indicate that metal-insulator transition 
in doped semiconductors occur at a critical density close to the Mott criterion [Ul ITT] . 
In the case of silicon and germanium doped with different dopants, the critical density 
varies as expected with the magnitude of Bohr radius with small, and may be signihcant, 
deviations [42] . 
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SYSTEM 

ao (nm) 

a^ (nm) 

nc2(cm ^) 

Si 

0.54 

1.5 

7.5 X 10^0 

Ge 

0.57 

3.7 

1.8 X 10^0 

SrTiOs 

0.39 

600 

1.5 X 10^^ 

PbTe 

0.65 

800 

4.6 X 10^6 


Table 2. The second critical doping estimated using, the lattice parameter and the 
effective Bohr radius of the four semiconductors. Above this threshold density, the 
local potential level is mostly set by contribution of global potential wells. 


Let us recall that when the equality expressed by Eq. 7 is fulhlled, metallicity 
emerges thanks to percolation. Yet many Bohr spheres do not touch each other (Fig.3b). 
For higher dopant concentration and/or much longer Bohr radius, we reach a regime 
with many Bohr spheres interpenetrating each other (Fig. 3C). Let us examine the 
consequences. 

4. Connection to the rongh Fermi seafloor 

Consider a lattice with a lattice parameter of ao, hosting a random distribution of wells. 
At an occupied site i, the local potential can be expressed as: 

r<ao^V, = Vo r > ao ^ V, = ( 8 ) 

r 

Here, Vq quantihes the maximum depth of the potential well, which can be 
assimilated to an isolated impurity level. Now, when there is a density n of such wells, 
the local potential level is set by the contribution of all occupied sites, near or far. The 
average potential level would be: 

< Vioc >= Voaon [ Airr^dr = AirVoaonr^p (9) 

Jo r 

The main new idea is that there can be threshold dopant density beyond which 
< Vioc > exceeds Vq. When this happens, the global contribution of potential wells 
weighs more than a single contribution event at the dopant site. In this regime, the 
global contribution shapes the local potential landscape. Such a situation occurs when 
the concentration attains a threshold value. Let us call this threshold nc 2 . It is given 
by expression: 

a*saon%^ = (- )V3 (10) 

71 

Note that this is expression is specific to three dimensions. A similar expression 
can be obtained for two dimensions by replacing Aur'^dr with 27 rrdr in the integral of 
Eq. 9. 

Table 2 gives the magnitude of nc 2 in the four semiconductors under discussion 
assuming that ao is the actual lattice parameter of the semiconductor. As seen in 
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the table, in both PbTe and SrTiOs, the magnitude of nc 2 is much lower than in 
silicon and germanium. It is tempting to attribute the striking difference in the density 
dependence of mobility seen in Fig. 2 to the fact that in contrast to the two elemental 
semiconductors, the long Bohr radius in PbTe and SrTiOa allows them to easily attain a 
regime were the local potential landscape is traced by the contribution of many potential 
wells. 

In a metal, where electrons can freely move, the chemical potential is everywhere 
the same but the local depth of the Fermi sea is not flat, because the environment 
created by dopants is random and not periodic (Fig. 3e). This real-space picture has 
its counterpart in the momentum space. The mobile electrons have a wave-vector and 
occupy a band originating from orbitals provided by atoms in the stoichiometric host. 
Therefore, there is a Fermi surface with a well-dehned locus in the Brillouin zone (F- 
point in the case of SrTiOs and L-point in the case of PbTe). The roughness of the 
seafloor reflects the stochastic distribution of dopants, which have generated metallicity. 
The distribution is assumed to be a Poisson distribution [ID], an assumption backed by 
experimental evidence |43]. 

In order to hnd the simplest model to describe the electron mobility in such a Fermi 
sea, we begin by recalling the simple expression linking the mobility of electrons to their 
mean-free-path, i: 

_e i 
^ hkp 

Now, what would be the intrinsic and elastic mean-free-path of electrons in such a 
context? This question can be replaced by another one: what is the average distance 
between two scattering centers in such a Fermi sea? A scattering center is where 
the (statistically unavoidable) deviation from the average density of dopants creates 
an island or a trench in the Fermi sea. Let us consider the mean-free-path divided 
by the interatomic distance — which is a dimensionless number. How many atomic 
sites shall one travel to hnd a place where the global contribution of dopants is such 
that there is a signihcant deviation from the average depth of the Fermi sea. In the 
simplest conceivable model, this number scales with two other dimensionless numbers, 
the interdopant distance and the rate of variation of a single potential well. Therefore: 



i _ n rxF 
— — o- 


( 12 ) 


In other words, this distance is simply proportional to the distance between dopants, 
n“^/^ and to the screening length, itf- The latter amplihes the scope of the well created 
by each extrinsic atom. C is a dimensionless parameter expected to be of the order of 
unity. Injecting the relevant expressions for n“^/^ and ttf-, one hnds an expression for 
intrinsic mobility in a Fermi sea generated from a random distribution of dopants. 


h-RD 


C e 
y/l^ao h 




(13) 
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According to this expression, mobility decreases with increasing carrier density in 
a given system. In different systems at the same carrier density, the longer the Bohr 
radius, the larger the expected mobility. Note that we are focusing on the mobility 
of isotropic three-dimensional crystals when their carrier density is larger than nc 2 as 
defined by Eq. 10. 

5. Comparison with experimental data 

Fig. 5 presents the mobility data reported in oxygen-deficient and Nb-doped SrTiOa by 
four different groups [331 IMl 1351IT^ . No disagreement is visible between the data reported 
by different groups. For comparison, we have also included the data reported on La- 
doped SrTiOajH], which presents a mobility slightly larger than Nb-doped samples in 
the dilute limit. As seen in figure, the expression for yiR£, gives a reasonable account 
of the data. The adjustable parameter C is roughly three times lower in the oxygen- 
deficient case. We will come back to this. 

Fig. 6 presents a similar treatment of the available data on PbTe[32l ITT] and 
BiTlSeS, which was the subject of a recent study. Novak et a/. [15] quantified mobility 
as a function of concentration in a system with a vanishing gap tuned half-way between 
BiTlSe 2 and BiTlS 2 , two band insulators (one trivial and one topological) with gaps 
of comparable amplitudes [IB]. The light cyclotron mass(m*=0.14me) combined to a 
relatively large dielectric constant(e ~ 80eo) yields an effective Bohr radius of a^ ~ 30 
nm. As seen in Fig. 6, the expression for /i^^) yielded by Eq.l3 seems a satisfactory 
description of the experimental data. Note the amplitude of parameter C, which is large 
compared to SrTiOs. 

Note that in this crude treatment, we have neglected the fact that neither the 
effective mass nor the dielectric constant do not stay constant as the carrier density 
is modified by several orders of magnitude. The apparent success of the fit indicates 
that the principal ingredients are present in the above-mentioned picture. Note that a 
blind linear fit of the logarithm of mobility vs. logarithm of density yields an exponent, 
which is somewhat lower than -5/6) (See table 3). Note that in all cases a downward 
deviation from the power-law behavior is visible at the lowest concentration. This is 
not surprising as the model is expect to brak down with the approach of nc 2 . 

What physical properties set the dimensionless fitting parameter C? To answer 
this question, let us consider two opposing effects of the change in carrier density on 
scattering rate. We have neglected them up to this point. The first is set by the 
comparative depth of the Fermi sea and the scattering potential. As the sea becomes 
deeper, the deviation from average concentration required for an island to emerge 
becomes harder to satisfy. Moreover, the deeper the potential well, the easier it is 
to satisfy this criterion. This effect is proportional to Vo/E^?. Second, electrons with 
a short wave-length are easier to scatter than those with longer ones. The scattering 
cross section scales with l/k|,. A fortunate accident for the simple model discussed here 
is that the density dependence of these two tendencies cancel each other out. More 
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n(cm'^) 

Figure 4. Mobility as a function of carrier density in Nb-doped (top) and reduced 
(middle) SrTiOa. The experimental data are reported in references [331 [Ml[SSI [Ill ?]• 
Two data points represent the data reported for La-doped SrTiOa)?]. Solid lines 
represent the expected variation of urd according to Eq. 13. In the bottom panel 
they are compared together. 
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Figure 5. Mobility as a function of carrier density in PbTe (top)|351[T3 and TlBiSSe 
(bottom) [3S] ■ Solid lines show the expected variation of yuD according to Eq. 13. 


quantitatively, they lead to the following expression for C: 
1 Ef 


c = 


(14) 


<^okp Vq 2m*alVo 

Thus, light mass and shallow potential enhance C and, consequently, the absolute 
value of mobility. Knowing the effective mass and the lattice parameter, Eq. 14 permits 
one to extract Vq from the magnitude of the fitting parameter, C. The results are given in 
Table 3. As seen in the table, as one would expect, the depth of potential well is much 
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SYSTEM 

a 

C 

ao (nm) 

m*(me) 

Vo (eV) 

E» (eV) 

SrTii,,Nb,03 

-0.73 

0.87 

0.39 

1.8 

0.16 

3.2 HZ] 

SrTiOs-^ 

-0.72 

0.25 

0.39 

1.8 

0.56 

3.2 HU 

PbTe (self-doped) 

-0.73 

15 

0.65 

0.07 

0.1 

0.2 Hi 

BiTlSeS (self-doped) 

-0.65 

4.7 

0.45 

0.14 

0.06 

<0.04 Hi 


Table 3. Parameters of the four systems under scrutiny, a represents the exponent 
found by htting the data to y, (x n°‘. C is the parameter extracted by htting the data 
to Eq. 13. Vo is extracted from C, the effective mass, m*, the lattice parameter ao 
and Eq. 14. Eg is the the band gap. 


lower in the small-gap semiconductors. The difference between oxygen-deficient and 
Nb-doped SrTiOa is also understandable. Substituting Ti with Nb introduces a single 
electron at the center of an octahedron. Oxygen vacancy, on the other hand, introduces 
two electrons and distorts two adjacent octahedra, digging a deeper potential well. 

6. Summary 

This paper puts under scrutiny an experimental fact already reported but rarely 
analyzed. In metallic semiconductors with a long Bohr radius, mobility is large and 
continuously decreases with increasing carrier density. Both these features are absent 
in silicon or germanium. According to the argument proposed here, this happens 
because in these metallic semiconductors the global distribution of dopants shapes the 
local seafloor. For this to happen, a threshold concentration of dopants is required. 
This second critical concentration becomes attainable only when the Bohr radius is 
long enough. With a set of most unsophisticated assumptions, one can formulate 
an expression for the intrinsic mobility of a metallic semiconductor with random 
distribution of dopants. In spite of neglecting subtle details, the expression is in good 
agreement with the data on four different semiconductors. 

This work is supported by Agence Nationale de la Recherche through the 
SUPERFIELD project. 
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